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Abstract 

A new series for the Madelung constant M is derived on the basis of a representation given by 

R. Crandall (1999). We are able to write M = C + S, where S" is a rapidly convergent series, and 

the constant C is fundamental: 

1 ln2 47r 1 T (i) T (|) 

C = + + )3' ^ -1.747564594 .... 

8 4^ 3 2^/2 7r3/2V2 

The remarkable result is that even if the S term be discarded, this constant C gives 10 good 

decimal places of M. This result advances the state of the art in the discovery of what Crandall 

has termed "close calls" to an exact Madelung evaluation. We indicate related identities and how 

this fundamental 10-digit accuracy might be further enhanced. 
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I. INTRODUCTION 

The Madelung constant has fascinated mathematicians and physicists ahke over the past 
century. It is physically related to the electrostatic interaction of a sodium ion with all other 
ions in a perfect NaCl crystal. The Madelung constant, M, can be written as 



where the sum over integers m, n and p runs from — oo to +00. A prime over the summation 
sign indicates that the term corresponding to the zero vector (m, n,p) = (0, 0, 0) be avoided. 
The sum in eq. fll.l|) is only conditionally convergent. One way to proceed is to interpret 
eq. ()l.l|) as a sum over expanding cubes (Borwein et. al 1985). Another approach is to 
interpret M as the analytic continuation of a sum — with (m^ + n"^ +p^)~* appearing with 
3f?(s) > 3/2 — to a definite value at s = 1/2. 

It has been an aim of many researchers to recast such lattice sums in terms of functions 
that decay exponentially fast for large r = a/ m'^ + ra^ + p"^. Equation p.l|) can be evaluated 
using Ewald method (Ewald 1920). The Ewald expansion certainly enjoys rapid decay, but 
error functions need be evaluated and this can be computationally problematic. Important 
work on Madelung constant and lattice sums has been carried out by Hautot (1975), Zucker 
(1976), Glasser and Zucker (1980), Borwein et. al (1985) and Crandall and Buhler (1987). 
Many of these historical treatments involve rapidly decaying sums of elementary functions. A 
fascinating account of Madelung constant, together with some new representations (such as a 
finite-domain integral representation), has been given by Crandall (1999). The present work 
uses that reference essentially as a starting point in a new quest to accelerate convergence. 

II. MADELUNG CONSTANT 

We start with the identity 




(1.1) 





2 



This identity can be found in Sperb (1996). A special case of the above identity is also found 
in Glasser and Zucker (1980). An easy proof may of the identity may be obtained with an 
application of the Possion summation theorem to the integral (Gradshteyn and Ryzhik 1965) 

We note an important aspect of the identity in eq. ^Z.\\ . With the help of this identity 
we can relate a summation in 2z/ + 1 dimensional space to a summation in 2v dimensional 
space. In fact taking z/ = one can easily obtain a generalization of Lekner's (1998) work 
on Coulomb sums for a triclinic cell. 

Putting x = 0, z/=l/2, r = a/ rn^ + rfi + and summing over integers m^n and p we 
obtain 

S= E E ; 2 2 2 -2 (2-3) 

k=—oo m,n,p 

E ; + 47r y ^ ' 



TT 



+ n2 +p2 ^^^^^ + + p2 exp (^27r ^m^Tn^Tp^ j - 1 
where we have used 

^i/2(x) = y^exp(-x). (2.4) 
Now, the second sum on the Ihs in eq. ()2.3|) can be obtained in a closed form, 



5 = hm V V ^ (2.5) 

a^O ^ ^ /c2 + ^2 _^ ^2 _^ p2 _^ Q,2 

fc=— oo m,n,p 

= ,i,„(v i-^"^"^' 

\k,m,n,p k 

^-^ k"^ + m? + + »2 \ ^-^ k"^ + 

k,m,n,p \ k 



The first sum on the rhs of eq. p.5j) may be evaluated with help of results given in Zucker 
(1984): 

> -7^ — T ^ o= n2, 2.6 

k,m,n,p 

while the second part is simply 



Thus we obtain 



M = — 



ln2 TT 
2 ~ ~ ~ 3 



(2.8) 



\m+n+p 



a/tti^ + + exp ( 27r^m^ + + ) — 1 



Equation ()2.8p is one of the important results of this paper. We will now derive an alternate 
expression for M again starting from the basic identity in eq. ()2.1|) . This time we substitute 
X = 1/2, which summing over m, n and p as before leads to 



-r 



T y — 



TT 



E 



m,n,p 



(2.9) 



-1 



,m+n+p 



y/m? + n'^ +p'^ exp y27c^/m?^i^r}?~^i^j + 1 
Again the rhs of eq. ()2.9|) can be obtained in a closed form. We start with 



fc=— oo m,n 



E E 

■oo m. 

E 



\ m+n+p 



^ -j^-^ m+n+p 



(2,10) 



- E 



1 



k,m,n'i 



The first sum on the rhs is again be obtained from the results contained in Zucker (1984): 

/ m+n+p 



E 



(2.11) 



The second sum on the rhs of Eq. (|2.1U| ) can be easily evaluated to be vr^. Thus we obtain 

' m+n+p 



M = V2-n + 2 



^p ^/m?~+~r^~+^ exp (2'IT^ym'^ + n'^ + p'^j + 1 



2.12) 



We can get a third and fast converging expression for M by taking average over M from 
eas. (0|l and (ET^ . We obtain 



M 



1 ln2 27r 1 

"4 ~ 27 ~ 71 



(2.13) 



2E 



m,n,p 



\ m+n+p 



+ 77,2+^2 



exp ( 47rA/m2~Ijr7^2~Ijr^ j _ I 
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This is not the whole story yet. We can now make use of a beautiful expression for M 
derived by Crandall (1999): 

r(|)r(i)v^ 



M = -2n + 



IT 



3/2 



(2.14) 



\ m+n+p 



^y m? + 'n? + exp i^TT^m? + + p^j + 1 
Taking the average over the two values of M from Eq. ()2.13|) and Eq. ()2.14|) we obtain 



M 



1 ln2 47r 1 T (|) T (|) 

8 47r 3 2v^ 7r3/2v^ 



(2.15) 



:-i 



, m+n+p 



a/ + + p2 exp ( Sn^Jm? + v? + — 1 



This representation of M is of the form C + S* of our Abstract. Even if the whole series part 
in Eq. ()2.15|) be ignored, we obtain 



1 ln2 47r^^^r(l)r(|) 



8 47r 3 2v^ 7t^/^V2 
1.747564594(7). 



(2.16) 



The value of M thus obtained is correct to an astounding 10 good decimals. 

The present method also leads to identities such as mentioned below. Eliminating M 
from Eq. (jTH)) and Eq. (jTT^ we obtain 

' . . 1 ln2 TT 1 

> — csch 27rv^ + + ] = 1 r=. 

^ Jm^ + n2 + p2 \ ^ 4 27r3y2 

rri,n,p V ' ' -f^ * 

Similarly combining eq. ()2.13j) and eq. ()2.14|l one obtains 



(2,17) 



-1 



.m+n+p 



m,n,p 



:csch (^4:71 a/ + 



^2 _|_ p2 



1 In 2 27r 

'8 ~ 4^"^ T 



(2.18) 



+ 



r(|)r(|) 



2%/2 7r3/2v/2 

Equations ()2.17|1 and ()2.18|1 can be represented in alternate form involving four dimensional 
sums. Another important identity that can be derived with the help of eq. ()2.H) and ()2.14|1 
is 

" MsJMsJ (2.19) 
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III. CONCLUSION 



We have given a new representation of the Madelung constant. The representation given 
here is very fast. Ignoring the series altogether, one still obtains Madelung constant correct 
up to ten decimal figures. Here I note that if one could obtain the sum 

/ -I \m+n+p 

tii (3 1) 

k:^n,p (fc + i)' + 16K + n2+p2)' 

in a closed from, then by combining that result with the one given here, one may obtain a 
representation of Madelung constant, where even dropping the series part, one would still 
obtain Madelung constant correct close to 20 decimal points. A probable approach may 
be through theta function identities involving q-series. However, we are not familiar with 
a systematic method of converting a product of theta function into a q-series, if one exist. 
It is hoped that the approach pointed out here may open up ways for researchers to look 
afresh on whether the Madelung constant can be expressed in a closed form. 

Acknowledgment s 

I am thankful to Prof. R. Crandall for helpful comments and mathematical advice. 



Electronic address: satst27@pitt.edu 



^ Borwein D, Borwein J and Taylor K 1985 J. Math. Phys. 16 1457-60. 

2 Ewald P P 1920 Die Berechnung optischer und elektrostatischer Gitterpotentiale, Ann. Phys. 
64 253-287 

^ Crandall R E 1999 Experimental Mathematics 8 367-379 
^ Hautot A 1975 J. Phys. A: Math Gen. 8 853-62 
5 Zucker I J 1976 J. Phys. A: Math Gen. 9 499-505 

^ Glasser M and Zucker I J 1980 Theoretical Ghemistry, Advances and Perspectives vol 5 (New 

York: Academic) pp 67-139 
^ Crandall R E and Buhler J P 1987 J. Phys. A: Math Gen. 20 5497-5510 

8 Sperb R 1998 Mol. Simulation 20 179-200 

9 Lekner J 1998 Mol. Simulation 20 356 



Zucker I J 1984 SIAM J. Math. Anal. 15:2 406-413 

Gradshteyn I S and Ryzhik I M 1965 Table of integrals series and products (Academic Press) 



7 



